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ABSTRACT 


The  results  of  four  digital  computer  calculations  based 
on  solutions  by  the  method  of  characteristics  of  the  equations 
for  stress  propagation  in  a  nonlinear  viscoelastic  material, 
with  a  plane  impulsive  source,  as  given  in  two  previous  reports, 
are  presented  in  graphical  form.  Many  of  the  observed  features 
of  the  variation  of  stress,  strain,  and  particle  velocity  are 
reproduced  by  the  theory,  including  an  increase  in  elastic 
constants,  as  compared  with  static  measurements,  and  the 
development  of  a  plateau  of  constant  residual  strain  near 
the  impact  face. 


LIST  OF  SYMBOLS 


cr  -  stress 

y  -  particle  velocity 

£  -  infinitesimal  strain  or  space  gradient  of  displacement 
£  -  Young ' s  modulus 
o^(r  -  dynamic  constants  of  material 
^  -  density 

C  -  velocity  of  initial  wave  front 
iz.  -  time 

?C  ~  Lagrangian  space  variable 
-  function  of  stress 
~  inverse  function  of  stress 
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INTRODUCTICW 


1  2 

The  two  previous  reports  ,  in  this  series  are  based  on 
Malvern'S  theory  j,  for  stress  propagation  in  strain-rate 
sensitive  materials.  They  are  concerned  with  the  loading  and 
unloading  parts  of  a  complete  stress  cycle  and  will  be  re¬ 
ferred  to  as  References  1  and  2.  Reference  1  deals  with  the  load¬ 
ing  process  and  extends  Malvern's  theory  to  include  finite- 
amplitude  strains  and  an  additional  dynamic  constant.  These 
modifications  provide  a  theoretical  basis  for  predicting  the 
residual  strain  distributions  and  the  increase  in  Young's 
modulus  that  are  observed  experimentally  from  the  dynamic 
testing  of  soils  and  of  some  metals.  Reference  2  is  primarily 
concerned  with  the  unloading  process  and  qualitatively  traces 
the  change  in  shape  of  a  rectangular  stress  pulse  applied  at 
X  =  0  as  it  propagates  into  the  material  j,  under  the  assumption 
of  a  linear  unloading  law.  The  theory  is  extended  to  include 
nonlinear  unloading  laws.  Numerical  procedures  are  given  for 
calculating  the  variation  of  stress,  strain,  and  particle 
velocity  during  the  complete  stress  cycle. 

This  report  presents  graphically  and  analyzes  the  results 
of  the  computations  performed  at  Ballistic  Research  Laboratories 
to  date.  One  set  of  calculations  is  based  on  static  stress- 
strain  data  for  the  Fort  Peck  Sand  as  determined  by  Whitman  , 


et  al^.  Two  sets  of  empirically  determined  dynamic  constants 

are  utilized  in  these  calculations,  one  based  on  a  value 

suggested  by  Parkin^;  the  other  is  empirically  determined  to 

compare  qualitatively  with  dynamic  tests  on  copper  reported 

7 

by  Kolsky  and  Oouch  ,  The  calculation  is  repeated  with  the 
second  set  of  dynamic  constants,  utilizing  a  different  form 
for  the  static  stress=strain  law. 

Another  computation,  based  on  physical  parameters  repre¬ 
sentative  of  dune  sand,  is  compared  to  numerical  results  from 
an  approximate  analytical  solution  derived  in  Ref.  1.  This 
comparison  shows  that  over  a  considerable  range  of  time  and 
distance,  the  approximate  analytical  solution  yields  results 
within  several  per  cent  of  the  solution  calculated  by  the 
digital  computer. 

Based  on  the  calculations  for  the  Fort  Peck  Sand,  the 
residual  strain  distribution  under  the  assumption  of  a  linear 
unloading  law  of  the  form  £-  j is  evaluated  for  various 
pulse  widths.  There  is  qualitative  agreement  between  these 
calculations  and  residual  strain  distributions  observed  in 
copper  when  the  rise  time-to-pulse  width  ratio  is  approximately 
the  same  between  theory  and  experiment . 
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FORT  PECK  SAND 


Physical  Constants 

The  static  stress=strain  response  of  the  Fort  Peck  Sand, 
experimentally  determined  by  Whitman,  et  al^,  is  shown  in 
Figure  1.  Young's  modulus  measured  as  the  initial  slope  of 
this  curve  is  (14.7)  10®  dynes/cm^,  and  the  density  is  re- 
ported  as  1.60  gm/cm^.  An  analytical  representation  of  the 
static  stress-strain  curve  is  obtained  by  fitting  an  empirical 

Q 

formula  suggested  by  Osgood  of  the  form  (Ref.  1,  eq.  2) 

£  =  f  ^  ^  (f)^ 

to  the  static  data.  The  constants  R  and  S  in  (1)  may  be 
determined  by  plotting  log  (  S.  -  cr/£-  )  vs.  log  (  cT/fi-  )  and  fitting 
a  straight  line  through  the  resulting  points.  The  slope  and 
intercept  of  this  straight  line  are  S  and  log  R,  respectively. 
The  constants  R  and  S  evaluated  in  this  manner  are  (6.86)  10^ 
and  2.22.  With  these  values,  (1)  adequately  fits  the  data 
over  the  range  of  interest  and  is  shown  superimposed  on  the 
experimentally  determined  stress-strain  curve  in  Figure  1. 

With  (1)  as  the  static  stress-strain  law,  the  basic 
constitutive  equation  defining  the  material  (Ref.  1,  Eq.  10) 
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takes  the  form 


i  =  i;[£- g-ccr^j  (2) 

If  this  is  interpreted  as  an  elastic-plastic  law  as  in 
Malvern's  work,  the  elastic  strain  rate  is  » 

while  the  simultaneously  occurring  plastic  strain  rate  is 

s't  =  £t-  £.=  =  i 

Thus^  the  plastic  strain  rate  is  assumed  proportional  to  the 
difference  between  the  instantaneous  strain  and  the  strain 
which  would  result  if  the  instantaneous  stress  had  been  stati¬ 
cally  applied.  In  Malvern’s  theory,  the  plastic  strain  rate 
was  a  function  of  the  difference  between  the  instantaneous 
stress  and  the  stress  the  material  would  support  if  the  instan 
taneous  strain  had  been  produced  statically. 

To  gain  some  insight  into  the  difference  between  these 
alternative  expressions  for  the  plastic  strain  rate,  calcu¬ 
lations  were  performed  using  identical  boundary  conditions 
and  physical  constants- -first  using  (1)  as  the  static  law  and, 
second,  solving  (1)  for  CT  so  that  (2)  was  proportional  to 

There  are  no  experimental  data  known  to  the  author  from 
which  the  dynamic  constants  a  and  b  may  be  directly  determined 
However,  it  is  known  from  theory  that  b/a  ^  E.  On  the  basis 


=  4- 


of  this  relation,  b/a  is  arbitrarily  chosen  as  (  14.9)  10®  dyne/cm^, 
which  yields  an  initial  wave  front  velocity  of  12,000  in/sec. 

One  set  of  dynamic  constants  is  found  by  assuming  l/a  =(6.5)10^- 
as  suggested  by  Parkin^  from  an  empirical  study  of  Whitman's 
data.  This  yields  b  =  3.32308.  The  calculation  using  these 
constants  and  (1)  as  the  static  law  is  referred  to  as  Case 
This  value  of  b  applies  to  a  hypothetical  material  that  ap¬ 
proaches  its  equilibrium  state  very  slowly,  as  is  illustrated 
by  the  fact  that  at  x  s  o,  almost  7  seconds  are  required  for 
the  material  to  reach  90%  of  the  ultimate  strain.  Since  such 
slow  flow  is  not  observed  experimentally  in  soils  and  metals, 
another  value  for  b  may  be  obtained  by  requiring  the  strain 
to  reach  90%  of  its  ultimate  value  in  0.005  second.  By  setting 
(9,(305^  =  .9  (^  (  )  and  solving  for  b  (Ref.  1,  Eq.  42), 

one  obtains 

=(2.38969)10“3 

By  requiring  the  initial  wave-front  velocity  to  be  12,000  in/sec, 
one  obtains  a  =  G/jsC^  =(  1 . 10633)  10‘^.  The  calculation 
using  these  more  realistic  values  for  a  and  b  with  (1)  as  the 
static  law  is  referred  to  as  Case  II .  The  calculation  with  the 
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dynamic  constants  for  Case  II  and  the  static  law  written  as 
is  referred  to  as  Case  III . 

Step  Response 

The  response  of  the  Fort  Peck  Sand  to  a  step  in  stress 
of  50  Ibs/in^  applied  at  x  =  0  is  numerically  evaluated  for 
the  above  two  sets  of  dynamic  constants  and  two  different  forms 
for  the  static  law.  These  solutions  are  compared  to  the 
centered  simple  wave  solution  obtained  by  neglecting  time-rate 
effects.  The  present  formulation  does  not  distinguish  between 
an  initial  tensile  or  compressive  stress.  A  tensile  stress 
is  assumed  positive  along  the  outward  normal  corresponding 
to  the  negative  x  direction.  The  resulting  particle  velocity 
is  parallel  to  the  stress,  but  is  negative  in  sign  since  the 
velocity  vector  is  determined  relative  to  the  coordinate 
system.  A  positive  velocity  would  correspond  to  a  compressive 
stress . 

Contours  of  constant  stress,  strain,  and  particle 
velocity  using  the  set  of  dynamic  constants  of  Case  I  (Figs.  2, 
3,  and  4)  show  the  initial,  slowly-decaying,  discontinuous 
elastic  wave-front,  followed  by  a  region  in  which  the  stress, 
strain,  and  particle  velocity  increase  toward  their  respective 
limits.  In  this  region,  lines  of  constant  amplitude  are  not 
straight;  however,  their  deviation  from  linearity  is  small 
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over  a  considerable  range  of  amplitude  and  strain  rate.  Because 
of  this,  the  resulting  contour  patterns  resemble  a  simple  wave 
originating  from  some  unknown  origin. 

The  calculation  is  repeated  using  the  second  set  of 
dynamic  constants  (Case  II),.  These  constants  essentially  con¬ 
tract  the  time  scale  compared  to  Case  I,  resulting  in  a  much 
closer  approximation  to  the  centered  simple  wave  solution,  and 
in  better  agreement  with  strain  records  obtained  from  a  soil 
shock  tube.  The  resulting  contours  of  stress,  strain,  and 
particle  velocity  (Figs.  5,  6,  and  7)  show  that  the  initial, 
discontinuous  wave-front  decays  quite  rapidly,  subsiding  to 
half  its  initial  amplitude  after  approximately  14  inches. 

As  in  Case  I,  lines  of  constant  amplitude  behind  the  initial 
wave  front  are  not  quite  straight.  These  limited  calculations 
suggest  that,  except  for  the  infinite  strain  rates  associated 
with  discontinuities,  phase  velocities  are  relatively  more 
sensitive  to  strain  amplitude  than  to  strain  rate. 

As  indicated  above,  the  Case  III  calculation  (Figs.  8, 

9,  10)  differs  from  the  Case  II  calculation  only  in  assuming 
the  plastic  strain  rate  proportional  to  the  stress  difference 
rather  than  the  strain  difference.  Comparison  of  the  contour 
patterns  for  these  two  cases  shows  that  the  material  approaches  equi¬ 
librium  much  more  rapidly  when  the  plastic  strain  rate  is  taken  pro- 
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portional  to  the  strain  difference  .  This  may  be  explained  by 
interpreting  the  characteristic  equations  for  a  nonlinear 
viscoelastic  material  (Ref.  1,  Eqs.  23,  24,  and  25)  as  a  pertur¬ 
bation  of  the  corresponding  plastic  equations  and  noting  that 
the  perturbing  term  (  ^(JT)-  )  is  in  general  numerically 

smaller  when  written  as  a  strain  difference  /  i •©•  >  ""  £.j_/ 

than  when  written  as  a  stress  difference  ^i.e.,  ( — /“ 
As  the  perturbing  term  approaches  zero,  the  solution  should 
approach  the  simple  wave  solution. 

For  purposes  of  comparison,  the  centered  simple  wave 
solution  based  on  the  static  stress-strain  law  (1)  is  shown 
in  Figure  11.  Since  it  is  based  on  the  rate-independent 
plastic  wave  theory,  it  may  be  directly  compared  to  the  calcu¬ 
lations  for  all  three  cases. 

Residual  Strain  Distributions 

q  Q 

A  commonly  accepted  procedure  ,  for  calculating  the 
residual  strain  distribution  is  to  assume  linear  unloading 
parallel  to  the  initial  slope  of  the  loading  cycle.  Under 
this  assumption,  the  residual  strain  distribution  resulting 
from  a  rectangular  stress  pulse,  t^  seconds  in  duration,  is 
(Ref.  2,  eq.  13) 
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with  the  loading  cycle  calculated  from  the  rate-independent 
theory,  the  residual  strain  produced  by  a  rectangular  stress 
pulse  of  variable  duration  shows  the  familiar  plateau  of 
constant  strain  near  the  impact  boundary  x  =  0  (Fig.  12). 
Beyond  this  plateau,  the  residual  strain  rapidly  approaches 
zero. 

Residual  strain  distributions  for  Case  I,  determined 
by  applying  (4)  to  the  Fort  Peck  Sand,  are  shown  in  Figures  13 
to  15  for  a  rectangular  stress  pulse  of  duration  1,  4,  and  10 
seconds  respectively.  These  results  differ  markedly  from  those 
based  on  the  rate°-independent  theory  since  no  plateau  of 
constant  residual  strain  is  developed. 

The  residual  strain  distributions  based  on  the  Fort 
Peck  Sand  (Case  II)  are  shown  in  Figures  16  to  19  for  pulses 
of  10,  20,  30  and  40  milliseconds,  respectively.  A  plateau 
of  almost  constant  strain  develops  close  to  x  =  0  for  a  pulse 
duration  which  is  several  times  greater  than  the  rise  time. 

For  this  calculation,  the  rise  time  is  arbitrarily  chosen  as 
the  time  required  to  attain  90%  of  the  ultimate  strain^-  0.005 
second.  Thus  the  observed  plateau  of  constant  residual  strain 
is  predicted  in  viscoelastic  theory,  if  appropriate  values  of 
the  dynamic  constants  are  chosen  and  if  the  pulse  is  long 


enough  for  the  material  near  the  impact  face  to  approach 
equilibrium. 

Residual  strain  distributions  resulting  from  impact 
tests  on  copper  bars  have  been  reported  by  Kolsky  and  Douch 
(Fig,  20).  The  dynamic  constants  for  the  Case  II  calculation 
were  determined  empirically  in  an  attempt  to  reproduce  these 
residual  strain  distributions.  This  was  done  by  observing 
that  the  ratio  of  rise  time  to  pulse  duration  for  the  experi¬ 
mentally  recorded  pulses  (Fig.  21)  which  produced  these 
residual  strain  distributions  was  of  the  order  of  four  to  five 
The  residual  strain  distributions  calculated  from  Case  II, 
where  the  corresponding  ratio  is  of  this  order  of 
magnitude,  are  qualitatively  similar  in  form  (Figs.  18  and  19) 
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DUNE  SAND 


Physical  Constants 

Typical  values  for  dune  sands  as  reported  by  Heiland^® 
are  =  1.76  gm/cm^,  1.64  ^  10“ ®E  ^  3.69  dynes/cm^.  The 
ratio  b/a  is  chosen  as  (3.69)  10®  dynes/cm^  and  E  as 
(1.64)  10®  dynes/cm^.  The  constant  b  is  calculated  by  assuming 
Malvern's  value  l/a  =  10®  for  metals,  yielding  b  =  7.7.  This 
calculation  was  performed  to  evaluate  the  validity  and  useful¬ 
ness  of  an  approximate  solution  (Ref.  1,  pp.  23-28),  which 
utilizes  a  static  stress-strain  curve  of  the  form  (1)  with 
S  =  2.0  for  mathematical  simplicity.  Another  restriction 
on  the  approximate  solution  requires  that  the  constant  R  be 
determined  from  the  relation 

2:  -  J-  ^  2  /f  o;  (®) 

(r-  '  r  ^  -ji 

With  CJ^  =  50  Ibs/ft^  as  one  boundary  condition  and  the  pre¬ 
viously  determined  constants,  R  =(5.54)io“®. 

Step  Response 

Contours  of  constant  stress,  strain,  and  particle  velocity 
from  digital  computer  data  again  show  a  simple  wave  pattern 
from  an  unknown  origin  (Figs.  22,  23,  and  24). 

The  solutions  calculated  by  the  digital  computer  are 
feplotted  to  facilitate  comparison  with  data  from  the  approxi- 
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mate  solution.  Figures  25"28  show  values  for  stress  and  strain 
as  functions  of  the  digitized  "reduced  time"  n,  calculated  from 
these  two  methods  at  representative  intervals.  The  arrival  of 
the  initial  discontinuous  wave  front  is  represented  by  n  =  0, 
and  each  integer  thereafter  represents  a  tisK  interval  of  0.05 
second.  Thus  n  =  15  corresponds  to  a  pulse  width  of  0.75  second. 

One  of  the  basic  assumptions  inherent  in  this  approxi¬ 
mate  solution  restricts  its  region  of  validity  to  small  "reduced 
times";  i.e.j  small  n.  As  can  be  seen  from  Figures  25  -  28, 
the  two  solutions  are  identical  at  n  =  0  and  diverge  slowly 
with  increasing  n.  The  percentage  errors  for  stress  and  strain 
over  the  range  of  available  data  have  been  calculated  and  the 
contours  of  constant  percentage  plotted  (Figs.  2Q  and  30). 

These  figures  show  that  the  maximum  error  for  stress  over  the 
range  of  data  examined  is  slightly  over  0.3%.  The  maximum 
error  for  strain  over  this  same  range  of  data  is  slightly  over 
7%.  This  calculation  is  deliberately  extended  beyond  its 
expected  range  of  application  since  n  =  15  corresponds  to  a 
loading  pulse  considerably  longer  in  duration  than  an  actual 
explosion  or  dynamic  test.  These  approximate  solutions  should 
be  valid  within  a  few  per  cent  over  the  range  of  pulse  durations 
normally  encountered.  The  use  of  other  values  of  a  and  b 
should  not  greatly  alter  the  general  conclusions  outlined  above. 
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COMPARISON  WITH  PLASTIC  WAVES  IN  METALS 


Stress-Strain  Curves 

Kolsky  and  Douch  have  determined  the  dynamic  and  static 

stress-strain  curves  for  an  aluminum  alloy;,  pure  aluminum,  and 

copper  (Figs.  31,  32  and  33).  The  tests  on  an  aluminum  alloy 

(Fig.  31)  show  that  the  values  of  stress  and  strain  determined 

dynamically  scatter  evenly  about  the  static  curve.  This  alloy 

belongs  to  the  class  of  strain-rate  insensitive  materials  which 
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are  well-described  by  the  Von  Karman-Taylor  plastic  wave  theory. 

Similar  dynamic  tests  on  pure  aluminum  show  that  it  is 
strain-rate  sensitive  (Fig.  32).  However,  there  is  no  discerni¬ 
ble  difference  between  the  dynamic  and  static  Young's  modulus. 
This  may  be  interpreted  in  terms  of  the  nonlinear  viscoelastic 
theory  (Ref.  1)  as  requiring  only  one  dynamic  viscoelastic 
constant.  Thus  it  should  be  possible  to  describe  the  dynamic 
behavior  of  pure  aluminum  by  Malvern's  theory,  utilizing  only 
one  dynamic  constant. 

On  the  other  hand,  the  dynamic  and  static  stress-strain 
curves  for  pure  copper  show  a  small  but  significant  increase 
in  Young's  modulus  under  dynamic  test  conditions  (Fig.  33). 

This  increase  in  Young's  modulus  suggests  that  at  least  two 
dynamic  constants  are  required  to  describe  adequately  the 
dynamic  behavior  of  copper,  as  outlined  in  Reference  1. 


The  data  from  which  these  dynamic  stress^^strain  curves 
are  constructed  would  be  obtained  from  gages  located  at  various 
distances  from  the  impact  face.  Since,  for  strain*^ rate  sensi¬ 
tive  materials,  there  is  no  one-to-one  correspondence  between 
stress  and  strain,  the  stress-strain  curve  deduced  from  experi¬ 
mental  data  will  be  a  function  of  recording  position.  To  gain 
some  insight  into  the  importance  of  recording  position,  the 
stress-strain  curves  obtained  from  the  calculations  of  Cases 
I  and  II  for  the  Fort  Peck  Sand  are  shown  in  Figures  34  and 
35  for  various  distances  from  the  impact  face.  There  is  an 
apparent  increase  in  the  yield  point,  which  decreases  with 
distance  from  the  origin,  and  in  Young's  modulus.  Above  the 
apparent  yield  point,  the  dynamic  curve  roughly  parallels 
the  static  curve  giving  the  appearance  of  shifting  the  entire 
static  curve  upward.  The  calculation  for  Case  II  shows  that 
for  distances  greater  than  12  inches,  the  apparent  stress- 
strain  curve  varies  slowly  with  small  changes  in  recording 
position. 


SUMMARY 


The  theory  presented  previously  (Ref.  1)  for  stress 
propagation  in  nonlinear  viscoelastic  materials,  assuming 
infinitesimal  strain,  has  been  programmed  for  a  digital  com¬ 
puter.  The  resulting  calculations  have  been  used  to  predict 
the  mechanical  parameters  of  stress,  strain,  and  particle 
velocity  as  a  function  of  time  and  distance  from  an  impulsive 
plane  wave  source.  By  suitable  choice  of  the  dynamic  constants, 
the  solution  can  be  made  to  approximate  a  simple  wave.  Such 
observed  features  as  an  increase  in  Young's  modulus  and  a 
plateau  of  constant  residual  strain  at  the  impact  boundary 
are  reproduced  by  the  theory. 

Three  calculations  were  performed,  permitting  comparison 
of  the  effect  on  wave  shape  of  different  dynamic  constants. 

The  first  two  calculations  (Cases  I  and  II)  assumed  that  the 
plastic  strain-^rate  is  proportional  to  the  difference  between 
the  instantaneous  strain  and  the  strain  which  would  result  if 
the  instantaneous  stress  had  been  statically  applied.  The 
dynamic  constants  of  Case  I  were  altered  for  the  calculations 
of  Case  II,  to  give  predictions  in  closer  agreement  with 
experimental  data.  In  the  third  calculation  (Case  III)  the 
plastic  strain  rate  was  taken  proportional  to  the  difference 
between  the  instantaneous  stress  and  the  stress  required 
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to  produce  the  instantaneous  strain  statically.  These  calcu¬ 
lations  indicate  that  the  former  assumption  results  in  a  more 
rapid  approach  to  values  that  would  be  obtained  under  static 
conditions . 

The  results  of  an  approximate  analytical  solution  (Ref. 

1)  are  compared  to  the  numerical  solution  from  the  digital 
computer.  Although  dynamic  constants  chosen  for  this  calcu¬ 
lation  were  not  physically  meaningful ,  the  results  indicate 
that  the  approximate  solution,  within  its  limits  of  validity, 
yields  results  which  are  sufficiently  exact  over  the  range 
of  time  and  distance  of  interest  for  most  applications. 

From  the  work  of  Kolsky  and  Douch ,  metals  appear  to  fall 
into  three  classes,  depending  upon  the  number  dynamic  constants 
required  to  describe  their  behavior.  Materials  such  as  certain 
aluminum  alloys  are  not  strain-rate  sensitive  and  thus  no 
dynamic  constants  are  required  for  describing  their  behavior 
under  impact.  Pure  aluminum  requires  only  one  dynamic  constant, 
since  Young's  modulus  is  the  same  under  static  and  dynamic 
conditions.  On  the  other  hand,  pure  copper  appears  to  require 
two  dynamic  constants  to  account  for  the  increase  in  Young's 
modulus  observed  under  dynamic  conditions. 
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